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ABSTRACT 


The  problem  under  consideration  is  that  of  estimating  simultaneously  the 
differences  between  the  means  of  p  >  2  test  treatments  and  the  mean  of  a 
control  treatment  when  the  population  variances  of  all  p  +  1  treatments  are 
known.  Tables  are  given  which  permit  the  experimenter  to  find  the  minimal  total 
number  of  experimental  units,  and  the  optimal  allocation  of  these  units  among 
the  p  +  1  treatments,  in  order  to  make  one-sided  or  two-sided  confidence 
interval  estimates  of  the  differences  of  interest.  These  intervals  achieve  a 
specified  joint  confidence  coefficient  1-a  for  a  specified  allowance  or 
"yardstick"  associated  with  the  common  width  of  the  intervals.  The  computations 
for  these  tables  are  based  on  the  results  of  Bechhofer  (1969)  for  one-sided 
comparisons,  and  Bechhofer  and  Nocturne  (1972)  for  two-sided  comparisons. 

Key  Words  and  Phrases:  Multiple  comparisons  with  a  control,  Dunnett's  procedure 
optimal  allocation  of  observations,  one-sided  comparisons,  two-sided  comparisons 
joint  confidence  coefficient,  completely  randomized  design. 
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1.  INTRODUCTION 

The  problem  of  comparing  simultaneously  p  >.  2  test  treatments  with  a 
control  treatment  arises  frequently  in  applied  research.  Dunnett  (1955), 

(1964)  considered  this  problem  and  provided  constants  necessary  to  make  joint 
lOO(l-a)  percent  confidence  statements  (either  one-sided  or  two  sided) 
between  the  mean  of  each  of  the  test  treatments  and  the  mean  of  the  control 
treatment  when  the  common  variance  of  the  p  +  1  treatments  is  unknown. 

An  important  design  decision  in  this  problem  is  that  of  how  to  allocate 
the  experimental  units  optimally  among  the  test  treatments  and  the  control 
treatment  when  the  p  t  1  variances  are  known  and  possibly  unequal.  Bechhofer 
(1969)  (hereinafter  referred  to  as  Bl)  gave  a  solution  to  this  problem  for 
one-sided  comparisons;  this  solution  is  globally  optimal  if  the  variances  of 
the  p  test  treatments  are  equal,  and  optimal  in  a  restricted  sense  if  these 
variances  are  unequal.  Bechhofer  and  Nocturne  (1972)  (hereinafter  referred 
to  as  B2)  generalized  these  results  to  two-sided  comparisons.  Bechhofer 
and  Turnbull  (1971)  gave  a  globally  optimal  solution  to  this  problem  for 
one -sided  comparisons  when  the  p  test  variances  are  known  and  unequal. 

Only  small  illustrative  sets  of  tables  of  optimal  allocations 
(aJl  for  p  =  2)  were  given  in  Bl  and  B2.  In  the  present  paper  we  give 
an  extensive  set  of  tables  for  p  =  2(1)10  both  for  joint  one-sided  or  joint 
twc -sided  comparisons  based  on  the  formulae  given  in  Bl  and  B2.  (See 
Remark  2.2  for  the  case  p  =  1.)  For  such  comparisons  these  tables  can  be 
use  d  to  determine  the  smallest  total  number  of  observations  necessary  to 
guarantee  selected  confidence  coefficients  (0.75(0,05)0.95,  0.99)  for  given 
specified  allowance  or  "yardstick"  associated  with  the  common  "width"  of  the 
confidence  intervals;  the  tables  also  tell  how  to  allocate  these  observations 
optimally  to  the  p  +  1  treatments. 
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Remark  1.1:  The  present  paper  (and  all  of  the  aforementioned  paper's)  deals  with 
the  case  in  which  a  completely  randomized  design  is  to  be  employed.  However, 
many  practical  situations  may  require  the  blocking  of  experimental  units.  If 
the  block  size  is  large  enough  to  accomnodate  one  replication  of  all  of  the  test 
treatments  and  additional  control  treatments  as  well,  then  the  optimal  allocations 
in  the  present  paper  can  be  used.  If  the  blocks  have  a  common  size 
k  <  p  +  1,  i.e.,  if  the  p  +  1  treatments  are  to  be  compared  in  incomplete 
blocks,  then  entirely  new  considerations  are  required  to  determine  the  optimal 
incomplete  block  design.  This  problem  is  considered  in  Bechhofer  and  Tamhane 
(1981). 

In  Section  2  we  introduce  our  notation  and  pose  the  optimal  allocation 
problem  both  for  one-sided  and  two-sided  comparisons.  The  tables  along  with  an 
explanation  of  how  they  are  to  be  used  are  given  in  Sections  3  and  4.  Section 
5  quantifies  the  loss  incurred  if  equal  allocation  is  used  instead  of  the 
optimal  allocation.  The  formulae  used  in  the  computation  of  Tables  1  through  9, 
and  details  of  the  computations  are  given  in  the  Appendix. 


2.  NOTATION  AND  STATEMENT  OF  THE  OPTIMAL  ALLOCATION  PROBLEMS 

Let  the  treatments  be  indexed  by  0,1,..., p  with  0  denoting  the  control 
treatment  and  l,2,...,p  denoting  the  p  >  2  test  treatments.  We  assume  that 

the  observations  (j  =  1,2,...)  on  the  ith  treatment  are  normally  dis- 

2 

tributed  with  unknown  mean  p.  and  known  variance  cr.  (0  <  i  <  p),  and  that 

l  i  —  — 

all  observations  are  mutually  independent.  Based  on  >.  1  observations  on 
the  ith  treatment  (0  <  i  <  p)  it  is  desired  to  make  either 


(I)  A  100(l-a)  percent  joint  one-sided  confidence  statement  of  the  form 
{y0  ’  Pi  =  X0  “  Xi  +  d  ^  1  -  P>K 


(2.1) 
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or 


(II)  A  100(l-cO  percent  joint  two-sided  confidence  statement  of  the  form 


(xQ  -  x.  -  d  <  iiQ  -  -  x.  +  d  (1  <  i  <  p)}. 


(2.2) 


In  (2.1)  and  (2.2),  x.  is  the  observed  value  of  the  random  variable 
Ni  1 

X.  =  J  X. ./N.  (0  <  i  <  p),  and  d  >  0  is  a  specified  allowance, 
l  .L.  13  i  “  =  r 

:=i 

The  optimal  allocation  problem  is  that  of  finding  the  allocation  vector 

2  2  2 

(Nq,Ni . N  )  which  for  known  (0^,0^, ... ,0^)  and  specified  1-a  and  d, 

P 

minimizes  the  total  sample  size  N  =  £  N.  subject  to 

i=0  1 


P{pn  -  u.  <  X  -  X.  +  d  (1  <  i  <  p)}  >  1  - 


(2.3) 


for  one-sided  comparisons,  and 


P{XQ  -  X.  -  d  <  uQ  -  u.  <  XQ  -  X.  +  d  (1  <  i  <  p)}  >l-o 


(2.4) 


for  two-sided  comparisons.  For  both  cases  we  denote  the  optimal  allocation  by 

A  A  A  A  P  •• 

(N. ,N, , . . . ,N  )  and  the  smallest  total  sample  size  by  K  =  I  (the  partic- 

U  J.  D  .  _ 


i  =  0 


ular  case  under  consideration  being  clear  from  the  context). 


Remark  2.1:  This  same  optimal  allocation  (NQ ,N^ , . . . ,N^ )  maximizes  the  joint 

2  2  2 

confidence  coefficient  for  known  (aQ ,a^ , . . . ,o  )  and  specified  total  sample  size 


Continuous  approximations  to  the  probabilities  (2.3)  and  (2.4)  are 
obtained  in  B1  and  B2,  respectively,  by  letting 


P 

y.  =  N./  J  N.  (0  <  i  <  p),  (2.5) 

1  1  i=0  1 


and  regarding  the  y.  as  nonnegative  continuous  variables  satisfying 
P  1 

y  Y.  =  1.  The  solutions  given  in  B1  and  B2  give  optimal  allocations 
i=0  1 
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for  one-sided  and  two-sided  comparisons  under  the  restriction  that  the  stan¬ 
dard  errors  (1  <.  i  <_  p)  of  the  test  treatment  means  are  equal, 

2  2,  2  2  2 
i.e.,  o^/y^  =  °j/Yj  J*  j;l  £  i,j  p);  if  °i  =  CT2  =  ‘  ‘  =  °p’  then  the 

solutions  give  globally  optimal  allocations. 

Under  the  stated  restriction  and  using  the  continuous  approximation,  the 

probabilities  (2.3)  and  (2.4)  can  be  shown  for  given  p  to  depend  on 
2  2  2 

(a.,0. ,...,o  ),  N  and  d  only  through  y„* 

0  1  D  0 

X  =  d,¥/°0  (2.6) 

and 

B  =  (  l  )/o2q.  (2.7) 

i=l 

For  given  p  and  B,  and  specified  1-a  the  optimal  solutions  which  we  denote 

»  a 

by  (yq,A)  are  uniquely  determined.  The  simultaneous  equations  which  yield 
these  solutions  are  given  in  the  Appendix. 

Remark  2.2:  It  should  be  noted  that  for  p  =  1  the  globally  optimal  allocation 

for  one-sided  and  for  two-sided  comparisons  is  aQ/N0  =  cr^/N^.  Then 

N  =  {(a  +  a,  )z  /d}2  and  N.  =  Na./(o„  +  a, )  (i  =  0,1)  for  one-sided  com- 

0  1a  l  l  0  1 

pari  sons,  the  same  expressions  holding  for  two-sided  comparisons  with  z^jy 

replacing  z^;  here  z^  is  the  upper  a-point  of  the  standard  normal  distri¬ 
bution. 

3.  DESCRIPTION  OF  THE  TABLES 

Tables  1  through  9  give  values  of  (Yq»*)  both  for  one-sided  and 
two-sided  comparisons  for  p  =  2(1)10,  respectively.  The  tabulated  values  of 
Yq  are  correct  to  within  one  in  the  fourth  decimal  place  while  the  tabulated 
values  of  A  are  rounded  up  in  the  fourth  decimal  place  to  guarantee  a  joint 
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confidence  coefficient  >  1-a  for  the  tabulated  value  of  yq.  For  each 
value  of  p  the  tabulations  are  made  for  1-a  =  0.75(0.05)0.95,  0.99  and 

0  =  p/2,  p,  3p/2  and  2p.  From  (2.7)  we  see  that  the  tables  can  be  used  for 

2  2  2 

the  special  case  a,  =  . . .  =  a  =  a  (say)  (in  which  case  the  allocations  are 

1  p 

2  2 

globally  optimal)  when  0  =  coQ  for  c  =  1/2,  l,  3/2,  2.  In  particular,  the 

2  2  2 

8  =  p  column  can  be  used  for  the  special  case  o„  =  a  =  . . .  =  o  . 

0  1  p 

An  examination  of  the  tables  shows  that  for  fixed  p,  8  and  1-a,  we  have 
and  X  in  the  two-sided  case  always  greater  than  the  corresponding  yq  and 
X,  respectively,  in  the  one-sided  case.  In  both  cases  for  fixed  p  and  8  we 
have  yq  increasing  with  1-a  and  approaching  the  limit  1/(1  +  /b)  as  1-a 
approaches  unity  (and  hence  Yq/y^  -*•  ao^°i  ^ or  1  i  p)*  This  limiting 

result  has  been  proved  analytically  in  B1  and  B2  for  the  one-sided  and  two- 

.  2  2  2 
sid^d  cases,  respectively.  For  0=0  =  ...  =  a  this  gives  the  limiting 

0  1  p 

result  that  Yq/y^  -*■  *^p  (1  <  i  <  p)  which  leads  to  the  recommendation  of 

Dunnett  (1955),  pp.  1106-1107  and  (1964),  pp.  486-487. 

4.  USE  OF  THE  TABLES 

A  2 
The  tables  of  (y„5X)  are  to  be  used  as  follows:  p  and  a.  (0  <  i  <  p) 

0  1  = 

are  given  as  data  of  the  problem;  these  determine  8.  The  experimenter  specifies 
d,  i.e.,  his  allowance,  and  his  one-sided  or  two-sided  joint  confidence  coef¬ 
ficient:  1-a.  Then  p,  8,  and  his  one-sided  or  two-sided  1-a  determine 

(Yq»X).  The  smallest  total  sample  size  N  is  then  the  smallest  integer  > 

~  2  -  -  - 
(X  oQ/d)  .  The  optimal  allocations  are  given  by  NQ  =  yqN  (to  the  nearest 

-  '  "  2  2 

integer)  and  N^.  =  (N-NQ)a^/80Q  (to  the  nearest  integer)  for  (1  <*  i  <  p); 
these  approximate  integer  allocations  which  were  obtained  by  using  the  con¬ 
tinuous  approximations  will  be  very  close  to  the  exact  integer  all  orations 
if  N  is  large. 
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Table  1 


y 

Optimal  allocation  on  the  control  (y  ) 

and  associated  X  to  achieve  a  given  joint 
confidence  coefficient  1-a 


P  =  2 


1  -ft 

6 

P/2 

P 

3p/2 

2p 

0 . 7S 

1-sided 

0.4245 

2.0074 

0.3519 

2.4818 

0.3115 

2.8420 

0.2841 

3.1441 

2-sidcd 

0.4638 

2.9007 

0.3868 

3.5403 

0.3435 

4.0266 

0.3141 

4.4348 

0.80 

1 

1-sided 

0.4417 

2.3226 

0.3666 

2.8556 

0.3246 

3.2608 

0.2961 

3.6009 

2-sided 

0.4691 

3.1474 

0.3908 

3.8352 

0.3468 

4.3588 

0.3169 

4.7986 

0.85 

i  1-sided 

0.4559 

2.6885 

0.3788 

3.2904 

0.3354 

3.7486 

0.3060 

4.1334 

i 

2-sided; 

0.4745 

3.4441 

0.3950 

4.1897 

0.3503 

4.7581 

0.3198 

5.2357 

0.90 

1-sided 

0.4683 

3.1481 

0.3893 

3.8376 

0.3448 

4.3631 

0.3146 

4.8048 

2-sided 

0.4802 

3.8298 

0.3995 

4.6506 

0.3540 

5.277] 

0.3231 

5.8039 

1-sided 

0.4801 

3.8298 

0.3992 

4.6510 

0.3536 

5.2780 

0.3225 

5.8053 

0.')!. 

2- sided 

0.4866 

4.4228 

0.4046 

5.3600 

0.3583 

6.0763 

0.3268 

6.6791 

1- sided 

0.4916 

5.1150 

i 

0.4085 

6.1894 

0.3616 

7.0114 

0.3296 

7.7036 

0.99 

2-:. idl'd 

0.4939 

5.5882 

0.4103 

6.7571 

0.3630 

7.6520 

0.3309 

8.4058 

1/(1  ‘  /p> 
which  is 
lim  y()  for 

(l-a)->l 

0.5000 

0.4142 

0.3660 

0.3333 

y 


The  upper  entry  in  each  cell  of  the  body  of  the  table  is 
Yq  and  the  iower  entry  is  X 
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Table  2 


Opt  i  in  a  1 
and  .u;., 


y 


n  the  control  (y  ) 


al  loi  .it  ion 
ociut<  d  A  to  achieve  a  p.iven  ioint 
confidence  coeff  icient  1-a 


P  =  3 


1-a 

fJ 

r 

P/2 

P 

— 

3p/2 

2P 

0.75 

1-sided 

0.3567 

2.6508 

0.2920 

3.3507 

0.2568 

3.8812 

0.2333 

4.3271 

2-sided 

0.4032 

3.6100 

0.3322 

4.4951 

0.2931 

5.1675 

0.2669 

5.7314 

0 .  80 

1-sided 

0.3766 

2.9904 

0 . 3086 
3.7561 

0.2713 

4 . 3378 

0.2464 

4.8259 

2-r.  ided 

0.4099 

3.8746 

0.3371 
4.8148  ' 

1 

0.2971 

5.5299 

0.2702 

6.1301 

0.85 

1-sided 

0.39  36 
3.3832 

0.3227 

4.2266 

0.2837 

4.8682 

0.2576 

5.4067 

2-sided 

0.4168 

4.1925 

0.3423 

5.1989 

0.3013 

5.9654 

0.2738 

6.6092 

0.90 

1-sided 

0.4089 

3.8757 

0.3353 
4.8186  | 

0.2948 

5.5366 

0.2675 

6.1398 

2-sided 

0.4241 

4.6057 

0.3479 

5.6986 

0.3059 

6 . 5320 

0.2777 

7.2326 

0 .9f< 

1- sided 

0.4240 

4.6059 

0.3475 

5.6993 

0.3053 

6.5336 

0.2770 

7.2350 

2-sided 

0.4324 

5.2417 

i 

0.3543 

6.4690 

0.3112 

7.4064 

0.2822 

8.1950 

0.99 

1-sided 

1  0.4389 

5.9855 

0.3591 

7.3725 

0.3151 

8.4332 

0.2856 

9.3263 

2-sidod 

0.4420 

6.4953 

0.3613 

7.9931 

0.3169 

9.1395 

0.2871 

10.1049 

1/(1  +  /$) 
which  is 
lim  Yq  f°r 
(1-a  )-*■! 

0.4495 

_ i 

0.3660 

0.3204 

0.2899 

y 


The  upper  entry  in 
Yg  and  the  lower 


eri 


noh  cell  of  the  body  of  the  table 
try  is  A 


is 


Table  3 


Optimal  allocation  on  the  control  (y  ) 

and  associated  A  to  achieve  a  given  joint 
confidence  coefficient  1-a 
P  =  4 


1-a 

3 

!  p/2 

p 

3p/2 

2P 

0.'7'> 

1-sided 

0.3163 

3.1961 

0.2570 

4.0932 

0.2251 

4.7742 

0.2040 

5.3454 

2-sided 

0.3645 

4.2095 

0.2978 

5.3084 

0.2617 

6.1432 

0.2376 

6.8435 

C 

a; 

o 

1-sided 

0. 3366 
3.5555 

0.2735 

4.5245 

0.2395 

5.2607 

0.2169 

5.8785 

2-sided 

0.3716 

4.4886 

0.3031 

5.6486 

0.2659 

6.5307 

0.2411 

7.2711 

0.8!) 

1 -sided 

0.3543 

3.9706 

0.2879 

5.0246 

0.2521 

5.8263 

0.2281 

6.4994 

2-sided 

0.3791 

4.8244 

0.3085 

6.0574 

0.2703 

6.9963 

0.2448 

7.7850 

0.00 

1- sided 

0.3706 

4.4904 

0.3011 

5.6535 

0.2634 

6.5394 

0.2383 

7.2836 

2-sided 

0.3870 

5.2607 

0.3144 

6.5894 

0.2750 

7.6026 

0.2488 

8.4553 

0.95 

1-sided 

0.3868 

5.2609 

0.3139 

6.5904 

0.2743 

7.6047 

0.2480 

8.4575 

2-sided 

0.3960 

5.9326 

0.3211 

7.4107 

0.2805 

8.5396 

0.2535 

9.4894 

0.99 

1-sided 

0.4031 

5.7197 

0.3263 

8.3762 

0.2846 

9.6432 

0.2569 

10.7100 

2-sided 

0.4064 

7.2602 

0.3286 

9.0412 

0.2865 

10.4043 

0.2585 

11.5523 

1/(1  +  /e) 

which  is 
lim  yQ  for 

(1-a  HI 

0.4142 

0.3333 

0.2899 

0.2612 

^The 


upper  entry  in  each  cell  of  the  body  of  the  table 


Yq  and  the  lower'  entry  is  A 
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Table  4 

Optimal  allocation-^  on  the  control  (yq) 
and  associated  X  to  achieve  a  given  joint 
confidence  coefficient  1-a 
P  "  5 


1-a 

P/2 

P 

3p/2 

2p 

1 -sided 

0.2884 

3.6806 

0.2331 

4.7562 

0.2036 

5.5731 

0.1842 

6.2583 

0.75 

2-sided 

0.3367 

4.7403 

0.2735 

6.0322 

0.2395 

7.0138 

0.2170 

7.8374 

0.80 

1-sided 

0.3086 

4.0569 

0.2493 

5.2096 

0.2176 

6.0857 

0.1967 

6.8210 

2- sided 

0.3440 

5.0325 

0.2787 

6.3901 

0.2437 

7.4229 

0.2205 

8.2830 

1-sided 

0.3263 

4.4911 

0.2635 

5.7351 

0.2299 

6.6817 

0.2076 

7.4764 

0.85 

2-sided 

0.3516 

5.3835 

0.2842 

6.8203 

0.2481 

7.9148 

0.2241 

8.8342 

0.90 

1-sided 

0.3428 

5.0344 

0.2766 

6.3961 

0.2411 

7.4334 

0.2176 

8.3050 

2-:'>idod 

0.3597 

5.8397 

0.2901 

7.3805 

0.2528 

8.5556 

0.2281 

9 . 5434 

0 

1 -sided 

0.3594 

5.8401 

0.2896 

7.3819 

0.2521 

8.5582 

0.2273 

9.5475 

2- si ded 

0.3690 

6.5430 

0.2969 

8.2463 

0.2583 

9.5473 

0.2327 

10.6420 

0.3762 

7. 3681 

0.3021 

9.2662 

0.2624 

10.7182 

0.2361 

11.9406 

0.99 

2-sided 

0.3796 

7.9354 

0.3044 

9.9699 

0.2642 

11.5272 

0.2377 

12.8389 

i/r  +  /b) 

which  is 
lim  Yq  for 
(l-a  HI 

0.3874 

0.3090 

0.2675 

0.2403 

y 


The  upper  entry  in  each  cell  of  the  body  of  the  table  is 
Yq  and  the  lower  entry  is  X 
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Table  5 


Optimal  allocation-^/  on  the  control  (Yq) 
and  associated  A  to  achieve  a  given  joint 
confidence  coefficient  1-a 
p  =  6 


1-a 

8 

P/2 

P 

3p/2  | 

2p 

0.7V 

J.-sided 

0.2677 

4.1224 

0.2154 

5.3627 

0.1878 

6.3052 

0.1696 

7.0960 

2-sided 

0.3154 

5.2230 

0.2549 

6.6925 

0.2227 

7.8095 

0.2014 

8.7469 

0.80 

1-sided 

0.2874 

4.5136 

0.2311 

5.8357 

0.2013 

6.8410 

0.1816 

7.6848 

2-sided 

0.3227 

5.5266 

0.2601 

7.0663 

0.2268 

8.2380 

0.2048 

9.2218 

0.85 

1-sided 

0.3050 

4.9646 

0.2451 

6.3838 

0.2132 

7.4639 

0.1922 

8.3710 

2-sided 

0.3303 

5.8913 

0.2655 

7.5156 

0.2311 

8.7532 

0.2084 

9.7930 

0.90 

1-sided 

0.3215 

5.5289 

0.2580 

7.0734 

0.2242 

8.2502 

0.2020 

9.2391 

2-sided 

0.3384 

6.3655 

0.2714 

8.1010 

0.2358 

9.4249 

0.2123 

10.5380 

0 . 9 5 

.l-sided 

0.3381 

6.3659 

0.2708 

8.1027 

0.2350 

9.4282 

0.2114 

10.5429 

2-sided 

0.3477 

7.0968 

0.2781 

9.0067 

0.2411 

10.4658 

0.2168 

11.6936 

0.99 

^  1-sided 

0.3551 

7.9560 

0.2832 

10.0755 

0.2451 

11.6970 

0.2201 

13.0624 

2-sided 

0.3584 

8.5476 

0.2855 

10.8140 

0.2469 

12.5493 

0.2216 

14.0108 

1/(1  +  v'B) 
which  is 
lim  Y0  for 
(1-a  HI 

0.3660 

0.2899 

: 

0.2500 

0.2240 

y 


The  upper  entry  in  each  cell  of  the  body  of  the  table  is 

A  * 

Yq  and  the  lower  entry  is  X 


Table  6 


Optimal  allocation  on  the  control  (y 

and  associated  X  to  achieve  a  given  joint 
confidence  coefficient  1-a 
P  =  7 


1-a 

8 

P/2 

P 

3p/2 

2p 

0.7') 

1-sided 

0.2513 

4.5320 

0 . 2016 
5.9265 

0.1754 

6.9865 

0.1583 

7.8761 

2-sided 

0.2983 

5.6693 

0.2402 

7 . 3047 

0.2094 

8.5487 

0.1891 

9.5920 

0.80 

1-oi ded 

0.2706 

4.9366 

0.2168 

6.4170 

0.1885 

7.5431 

0.1698 

8.4886 

2-sided 

0.3055 

5.9833 

0.2452 

7.6929 

0.2134 

8.994.? 

0.1924 

10.0871 

0.8', 

1 -sided 

0.2880 

5.4027 

0.2304 

6.9855 

0.2000 

8.1905 

0.1801 

9.2026 

2-sided 

0.3130  1 

6.3604 

0.2506  ' 

8.1597 

0.2176 

9.5308 

0.1959 

10.6830 

0.90 

1-sided 

0.3042 

5.9860 

0.2431 

7.7010 

0.2107 

9.008L 

0.1895 

10.1056 

2-sided 

0.3211 

6.8509 

0.2563 

8.7680 

0.2221 

10.2307 

0.1997 

11.4607 

0.95 

1-sided 

0.3208 

6.8515 

0. 2557 
8.7700 

0.2213 

10.2345 

0.1988 

11.4663 

2-sided 

0.3304 

7.6079 

0.2629 

9.7101 

0.2273 

11.3165 

0.2040 

12.6683 

0.99 

l-:;ided 

0.3377 

8.4982 

0.2679 

10.8236 

0.2312 

12.6029 

0.2072 

14.1013 

2-sided 

0.3410 

9.1120 

0.2701 

11.5941 

0.2329 

13.4947 

0.2086 

15.0956 

1/(1  +  fft) 
which  is 
lim  Yq  for 

(l-a)-*l 

0.3483 

0.2743 

0.2358 

0.2109 

y 


The  upper'  entry  in  r.irh  cell  of  the  body  of  the  table  is 
Yq  and  the  lower  entry  is  X 


1 


Table  7 

y 

Optimal  allocation  on  the  control  (y^) 

and  associated  X  to  achieve  a  given  joint 
confidence  coefficient  1-a 

p  -  8 


1-a 

8 

l>/2 

P 

3p/2 

2p 

0.75 

1-sided 

0.2380 

4.9161 

0.1904 

6.4561 

0.1654 

7.6273 

0.1491 

8.6103 

2-sided 

0.2841 

6.0870 

0.2230 

7.8787 

0.1984 

9.2416 

0.1790 

10.3857 

1 

1 

0.80 

1-sided 

0.2569 

5.3328 

0.2052 

6.9627 

0.1780 

8.2030 

0.1603 

9.2444 

2-sided 

0.2912 

6.4105 

0.2330 

8.2802 

0.2023 

9.7038 

0.1822 

10.9000 

0.85 

1-sided 

0.2739  ! 

5.8129 

0.2184 

7.5499 

0.1893 

8.8727 

0.1702 

9.9840 

2-sided 

0.2987 

6.7991 

0.2382 

8.7630 

0.2064 

10.2601 

0.1856 

11.5182 

0.00 

1-sided 

0.2899 

6.4135 

0.2308 

8.2892 

0.1997 

9.7192 

0.1793 

10.9211 

2-sided 

0.3066 

7.3046 

0.2438 

9.3926 

0.2108 

10.9860 

0.1892 

12.3260 

0.8  5 

l-nidod 

0.3063 

7.3053 

0.2432 

9.3949 

0.2100 

10.9904 

0.1883 

12.3325 

2- si ded 

0.3158 

8.0853 

0.2502 

10.3684 

0.2158 

12.1133 

0.1934 

13.5818 

0.99 

1- sided 

0.3230 

9.0045 

0.2551 

11.5233 

0.2196 

13.4509 

0.1965 

15.0743 

2-sided 

0.3263 

9.6389 

0.2573 

12.3235 

0.2213 

14.3794 

0.1979 

16.1112 

1/(1  *  */B) 
which  is 
lim  y0  for 

(1-a  HI 

0.3333 

0.2612 

0.2240 

0.2000 

y 


The  upper  entry  in  each  cell  of  the  body  of  the  table  is 

A  * 

Yq  and  the  lower  entry  is  X 


Table  8 


Optimal  allocation 


y 


on  the  control 


<v 


and  associated  X  to  achieve  a  given  joint 
confidence  coefficient  l-o 


p  =  9 


]-U 

1  1  1 

v> 

P 

3p/2 

2P 

0.  v 

1  -  sided 

0.2269 

6.2794 

0.1811 

6.9578 

0.1572 

8.2346 

0.141  6 
9.3066 

id-  d 

0.2721 

6.4812 

0.2177 

8.4214 

0.1892 

9.8977 

0.1705 

11.1372 

0  .  h » > 

m 

0.2453 

5.7073 

0.1954 

7.4792 

0.1693 

8.8280 

0.1523 

9.9608 

m 

0.2791 

6.8135 

0.2226 

8.8352 

0. 1930 
10.3750 

0.1736 

11.6686 

0.66 

1- sided 

0. 2620 
6.2002 

0.2083 

8.0837 

0.1802 

9.5186 

0.1619 

10.7242 

2-sided 

0.2864 

7.2128 

0 . 2277 
9.3331 

0.196't 

10.9496 

0.1769 

12.3084 

0  .MO 

1 -sided 

0.2777 

6.8170 

0.2204 

8.8452 

0.1903 

10.3919 

0.1708 

11.6921 

2-sidod 

0.2942 

7.7324 

0.2331 

9.9824 

0.2012 

11  .”000 

0.1804 

13. 1444 

0.9  6 

m 

0.2939 

7.7331 

0.2325 

9.9851 

0.2004 

11.7048 

0.1796 

1  3 .  IMS 

■ 

0.3033 

8.5352 

0.2394 

10.9897 

0.206: 

12.866a 

m 

0.99 

1- sided 

0.3104 

9.4814 

0.2441 

12.1834 

0 . 209  1 

14.2615 

0.1874 

16.9932 

0.3136 
10.1 351 

0.2463 

13.0113 

0.211  < 
16.214.' 

0.1867 

17.0649 

1/(1  t  /fl) 
which  is 
lim  yQ  for 

(l-o) >1 

0 . 3204 

0.2500 

0.2139 

0.1907 

V 


The  upper  entry  in  each  cell  of  the  body  of  the  table  is 
y  and  the  lower  entry  is  X 
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Table  9 

y 

Optimal  allocation  on  the  control  (y.) 

*  u 

and  associated  X  to  achieve  a  given  ioint 
confidence  coefficient  l-o 

p  =  10 


!  1 

\  J 


i 


\ 


1-a 

9  i 

P/2 

P 

— 

3p/2 

2p 

0 . 74 

l-sided 

0.2174 

5.6252 

0.1731 

7.4359 

0.1501 

8.8138 

0.1351 

9.9708  ’ 

2-sided 

0.2618 

6.8558 

0.2089 

8.9378 

0.1813 

10.5225 

0.1632 

11.8532 

0.80 

1-sided 

0.2354 

6.0635 

0.1871 

7.9712 

0.1619 

9.4237 

0.1455 

10.6438 

2-sided 

0.2686 

7.1964 

0.2136 

9. 3633 

0.1849 

11.0140 

0.1662  I 

12.4009 

O.i', 

1- sided 

0.2517 

6.5684 

0.1996 

8.5919 

0.1725 

10.1338 

0.1548  i| 

11.4295  1 

2- sided 

0.2758 

7.6058 

0.2186 

9.8752 

0.1888 

11.6058 

0.1694 

13.0607  . 

0.(0 

1 -sided 

0.2672 

7.2003 

0.2114 

9.3742 

0.1823 

11.0323 

0.1634 

12.4263  ^ 

2-sided 

0.2835 

8.1385 

0.2239 

10.5431 

0. 1930 
12.3790 

"  1, 

0.1728  > 

13.9232  I 

0.4  4 

1 -sided 

0.2831 

8. 1393 

0.2233 

10.5461 

0.1921 

12.3844 

0.1719 

13.9310 

2-sided 

0 . 2924 
8.9621 

0.2300 

11.5801 

0.1977 

13.5916 

0.1767 

15.2663 

0.99 

1 -sided 

0.2993 

9.9338 

0 . 2346 
12.8102 

0.2011 

15.0121 

0 . 1  79  5 

16.8667 

2-sided 

0.3025 

10.6057 

0.2367 

13.6643 

0.2027 

16.0071 

0.1808 

17.9810 

1/(1  ♦  /ft) 
which  is 
lim  yQ  for 

(!-•)►! 

0.3090 

0.2403 

0.2052 

0.1827 

t 


y 


The  upper  entry  in  each 


cell  of  the  body  of  the  table  is 


yQ  and  the  lower  entry  is  X 


2 


1  5 


Numerical  examples:  suppose  tnat  p  =  d  ana  a^  =  l  (0  =<  i  ^  3);  then 
8  =  p  =  3.  If  one-sided  intervals  are  desired  with  d  =  0.5  and 

1-a  =  0.95  then  from  Table  2  we  find  that  y^  =  0.3475,  X  =  5.6993.  Hence 

N  =  [{(5.6993)1/0. 5}2]  =  [129.9]  =  130  and  N  =  46,  N  =  N  =  N  =  28. 

0  1  t  o 

2  2 

For  the  same  specification  with  =  1/2,  =  1  (1  <_  i  <*  3)  we  have 

*  A  A  1  * 

8  =  2p  =  6  and  hence  yQ  =  0.2770,  X  =  7.2350;  thus  N  =  105,  =  30,  = 

=  25.  These  calculations  give  an  indication  of  the  sensitivity  of  the 

2 

allocations  and  sample  sizes  to  rather  large  changes  in  oQ  and  the  ratio 

2  >  2  2  2 
o./a()  (1  <  i  <  3)  when  =  a 

2  2 

If  the  experimenter  is  prepared  to  assume  that  o^=o  (0  ^  i  <  p) 

2 

where  the  actual  value  of  a  is  unknown,  and  he  believes  that 
2  2  2  2  2 

a  <  a  <  o  where  o  <  a  are  known,  then  this  information  can  be  used  in 
L  U  L  U 

2  2 

designing  the  experiment,  e.g.,  acting  as  if  a  =  leads  to  a  conservative 

choice  of  N.  However,  after  the  experiment  has  been  conducted,  when  the 

.  .  2 
results  are  being  summarized,  the  common  unknown  a  should  be  estimated 

2 

using  the  pooled  data,  the  estimate  s  being  based  on  v  =  N  -  (p+1)  d.f. 
The  estimate  should  then  be  used  with  Dunnett's  (1955)  formulae  for  joint  con¬ 
fidence  statements  (analogous  to  (I)  and  (II)  of  (2.1)  and  (2.2), 
respectively) : 


(III)  A  100(l-a)  percent  joint  one-sided  confidence  statement 

{Vi„  -  li.  <  x  -  x.  +  t^  A/- r—  +  -4—  (1  <  i  <  p)) 

0  i-O  i  v,p,p  7  N  -  - 


(4.1) 


or 


(IV)  A  100(l-a)  percent  joint  two-sided  confidence  statement 
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{x  -  x.  -  t'^ 

0  i  v,p,p 


±-4-i,0 

»o  »i 


V. 

1 


Vi’ =Vf  fi  a 4 1 4 p)1' 


(a) 


(4.2) 


Here  t^”'  p  is  the  upper  a  equicoordinate  point  of  the  p-variate 

t-distribution  (p-variate  | t | -distribution)  with  d.f.  v  and  equal  corre- 

*  A  A  (Cl  ) 

lations  p  =  N  /(N  +N  );  tables  of  t  are  given  by  Krishnaiah  and 

101  v,p,p 

Armitage  (1966)  while  tables  of  t'^a  are  given  by  Hahn  and  Hendrickson 

v,p,p 


(1971). 

5.  LOSS  FROM  EQUAL  ALLOCATION  ON  ALL  p  +  1  TREATMENTS 

It  is  of  some  interest  to  determine  how  much  is  lost  in  terms  of  increased 

sample  size  if  the  experimenter  uses  equal  allocation  on  all  p  +  1  treatments, 

i.e.,  Nq  =  =  ...  =  Np  =  N/(p  +  1),  instead  of  the  optimal  allocation.  Tor 

this  purpose  we  define  the  relative  efficiency,  RE  =  N/N.  Here  N  and  N  arc 

the  total  sample  sizes  required  to  guarantee  the  same  joint  one-sided  (two-sided) 

confidence  coefficient  1-a  using  (2.3)  (using  (2.4))  for  given 
2  2  2 

(°0,Ol’  ’ * ' ’°p)  anc*  specified  d  when  the  optimal  allocation  and  equal  allo¬ 
cation,  respectively,  are  employed.  (Note  that  RE  <  1,  and  small  values  of 

RL  indicate  large  relative  savings  by  optimal  allocation  over  equal  allocation.) 

.  2  2 
We  now  determine  RE  for  the  important  special  case  a^=o  (0  <.  i  <  p). 

(Other  cases  can  be  determined  analogously.)  For  one-sided  comparisons  we 

have  (ignoring  the  integer  restrictions) 


Thus 


(6.1) 


(6.2) 


m 
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for  one-sided  comparisons.  Here  t^01^  ^/2  t*ie  uPPer  “-point  of  the 
distribution  of  the  maximum  of  p  equicorrelated  standard  normal  random 
variables  with  common  correlation  P  =  1/2:  the  values  of  t^a^  ,  have  been 
tabulated  for  selected  p  and  1-a  by  Gupta,  Nagel  and  Panchapakesan  (1973). 

(a  J 

For  two-sided  comparisons  »  p  1/2  *n  (5.1)  and  (5.2)  is  replaced  by 
where  t’^ot\/_  is  the  upper  a-point  of  the  distribution  of  the 
maximum  of  the  absolute  values  of  p  equicorrelated  standard  normal  random 

t  /  \  2 

variables  with  common  correlation  P  =  1/2:  the  values  of  ( t  „  ,  '1  have  been 

»,p,l/2' 

tabulated  for  selected  \>  and  1-a  by  Krishnaiah  and  Armitage  (1965). 


Entries  for  v  =  60  the  tables  of  Hahn  and  Hendrickson  (1971)  can  be  used 

.  .  1 (a ) 

as  conservative  appro\  mations  to  t  , 

“.P .1/2 

Gome  representative  HE  -  values  for  one-sided  comparisons  (similar 


results  would  be  obtained  for  two-sided  comparisons)  are  given  in  the  following 
table.  Values  of  RE  for  One-Sided  Comparisons 


1-a 

P 

2 

5 

10 

0.75 

0.9986 

0.9741 

0.9390 

0.9  5 

0,9818 

0.9101 

0.8433 

1 .  oy 

[  0.9759 

0.8890 

0.8121 

It  can  be  seen  that  the  relative  savings  using  the  optimal  allocation  increase 

with  p  and  1-a  as  would  be  expected;  also,  for  fixed  p  and  1-a  the 

~  2 
difference  N  -  N  is  directly  proportional  to  (o/d)  . 
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APPENDIX 


Formulae  for  optimal  allocation ,  and  details  of  computation 

Aj_l  F ormulae  for  optimal  allocation  for  one-sided  comparisons  (Reference  Bl) : 

Let  $(•)  and  $(•)  denote  the  standard  normal  distribution  and  density 
function,  respectively,  and  let  !  P)  denote  the  equicoordinate  k -variate 

standard  normal  distribution  function  with  common  correlation  p .  Then  the 

given  in  Tables  1  through  9  are  the  unique  solutions  of  the  follow¬ 
ing  simultaneous  equations  (A.l)  and  (A. 2): 


/  *p[(—  +  X)  (^)1/2]  d*(x)  =  l-o, 

/7  3 


(A.l) 


[(1-6)y2-2y+1]t$p_1[t | <1-y)/{2(1-y)+Y$}] 

-  %-W1- Y)g  .{t)4  |_  1-Y  i  .  0 

2(1-y)+Y0  ?  ^  P-2  LTt3(l-Y)+YB7  Ud-vl+va-1  °’ 


where 


(A. 2) 


T  =  XY  (Ci. 


$(i-y) 


Y+y6][2(1-y)+yB] 


V/2 

V  • 


A. 2  Formulae  for  optimal  allocation  for  two-sided  comparisons  (Reference  B2). 

Here  the  (yq»A)  given  in  Tables  1  through  9  are  the  unique  solutions 
of  the  following  simultaneous  equations  (A. 3)  and  (A. 4): 

00 

/  {#[(—■+  X)(-~)172]  -$  [(-2=  -  X)  (^S^')^72]}Pd|J>(x)  =  1-a,  (A. 3) 

✓y  p  »Y  p 


—  Hill  Jli  tlli  ii  ■  -» 


19 


a[(1-B)y  -2y+l]D 


,1/2 


( l-y+Y$) 


1/2 


-  o. 


(A. 4) 


[2(1-y)+y6] 


y3 


whei’e 


Dx  =  <J>p_1{-A1T1’  Tp  |  (1-y)/(2(1-y)+yB}  ) , 

d2  =  ‘J>p_2^"A2T2’  t2^1_y^/,^3^1_y^+y6^’ 

D3  =  <I>p_  2  ^  —  T  3  t3K1-Y)/{3(1-y)+y£}), 

1  /2 

Tp  =  t,  t2  =  t[(1-y+yB)/{3(1-y)+yB}]  ,  T3  =  h1  t2, 

A1  =  {2(1-y)+yB)/yB,  A2  =  {4(1-y)+yB)/y$,  and 

4>,  (a,b|p)  =  P{a  <  Z.  <_  b  (1  <.  i  <  k)}  where  the  Z.  are  standard 
normal  with  corr{Z^,Z^}  =  p  for  i  ^  j,  l^i,  j  <  k. 

A. 3  Details  of  computation 

The  IMSL  subroutine  ZSYSTM  was  used  to  solve  the  pairs  of  simultaneous 

equations  (A.l),  (A. 2)  and  (A. 3),  (A. 4).  The  stopping  criteria  used  in 

arriving  at  the  final  solutions  were  the  following:  (i)  the  difference 

between  the  left  hand  and  the  right  hand  sides  of  each  equation  is  less  than 
“6 

1  x  10  or  (ii)  in  two  successive  iterations  the  corresponding  trial  values 
of  Yq  and  A  do  not  differ  in  the  first  six  significant  digits. 

To  evaluate  a  quantity  of  the  form  4^(a.b|p)  (which  includes  <J>^(b|p) 
as  a  special  case  for  a  =  -  “)  the  following  iterated  integral  represen¬ 


n 


tation  (see  equation  (2)  of  Bechhofer  and  Tamhane  (1974))  was  used: 
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For  p  =  2  the  quantity  ^  (a,b|p)  reduces  to  4>(b)-$(a)  and 

♦  (a,b|p)  =  ♦  (a,bjp)  =  1.  Thus  the  evaluation  of  the  various  expressions 

p-2  u 

is  particularly  simple  for  p  =  2. 

To  evaluate  $(•)  the  formula  (26.2.17)  given  in  Abramowitz  and  Stegun 
(1964)  was  used;  this  formula  is  accurate  to  within  ±  7.5  x  10  .  The  Romberg 

quadrature  method  was  used  to  evaluate  the  various  integrals.  All  of  the 
calculations  were  done  on  a  CDC  6600  computer  at  Northwestern. 


The  tabulated  values  of  are  rounded  off  in  the  fourth  decimal  place 

while  the  values  of  X  are  rounded  up  in  the  fourth  decimal  place  (to  insure  a 
joint  confidence  coefficient  >.  1-a). 
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